Metallicity Distribution Function Research Articles

Strong outflows and inefficient star formation in the reionization-era ultrafaint dwarf galaxy Eridanus ii

ABSTRACT We present novel constraints on the underlying galaxy formation physics (e.g. mass-loading factor, star formation history, and metal retention) at z ≳ 7 for the low-mass (M* ∼ 105 M⊙) Local Group ultrafaint dwarf galaxy (UFD) Eridanus ii (Eri ii). Using a hierarchical Bayesian framework, we apply a one-zone chemical evolution model to Eri ii’s CaHK-based photometric metallicity distribution function (MDF; [Fe/H]) and find that the evolution of Eri ii is well characterized by a short, exponentially declining star formation history ($\tau _\text{SFH}=0.39\pm _{0.13}^{0.18}$ Gyr), a low star formation efficiency ($\tau _\text{SFE}=27.56\pm _{12.92}^{25.14}$ Gyr), and a large mass-loading factor ($\eta =194.53\pm _{42.67}^{33.37}$). Our results are consistent with Eri ii forming the majority of its stars before the end of reionization. The large mass-loading factor implies strong outflows in the early history of Eri ii and is in good agreement with theoretical predictions for the mass scaling of galactic winds. It also results in the ejection of >90 per cent of the metals produced in Eri ii. We make predictions for the distribution of [Mg/Fe]–[Fe/H] in Eri ii as well as the prevalence of ultra metal-poor stars, both of which can be tested by future chemical abundance measurements. Spectroscopic follow-up of the highest metallicity stars in Eri ii ([Fe/H] > −2) will greatly improve model constraints. Our new framework can readily be applied to all UFDs throughout the Local Group, providing new insights into the underlying physics governing the evolution of the faintest galaxies in the reionization era.

Stellar Metallicities and Gradients in the Isolated, Quenched Low-mass Galaxy Tucana

We measure the metallicities of 374 red giant branch (RGB) stars in the isolated, quenched dwarf galaxy Tucana using Hubble Space Telescope narrowband (F395N) calcium H and K imaging. Our sample is a factor of ∼7 larger than what is available from previous studies. Our main findings are as follows. (i) A global metallicity distribution function (MDF) with 〈[Fe/H]〉=−1.55−0.04+0.04 and σ[Fe/H]=0.54−0.03+0.03 . (ii) A metallicity gradient of −0.54 ± 0.07 dex Re−1 (−2.1 ± 0.3 dex kpc−1) over the extent of our imaging (∼2.5 R e ), which is steeper than literature measurements. Our finding is consistent with predicted gradients from the publicly available FIRE-2 simulations, in which bursty star formation creates stellar population gradients and dark matter cores. (iii) Tucana’s bifurcated RGB has distinct metallicities: a blue RGB with 〈[Fe/H]〉=−1.78−0.06+0.06 and σ[Fe/H]=0.44−0.06+0.07 and a red RGB with 〈[Fe/H]〉=−1.08−0.07+0.07 and σ[Fe/H]=0.42−0.06+0.06 . (iv) At fixed stellar mass, Tucana is more metal-rich than Milky Way satellites by ∼0.4 dex, but its blue RGB is chemically comparable to the satellites. Tucana’s MDF appears consistent with star-forming isolated dwarfs, though MDFs of the latter are not as well populated. (v) About 2% of Tucana’s stars have [Fe/H] < −3% and 20% have [Fe/H] > −1. We provide a catalog for community spectroscopic follow-up.

Visit Nearby Halo Substructures Using LAMOST DR9 MRS Data

Based on the LAMOST DR9 Medium Resolution Spectra (MRS) catalog, combined with astrometries from Gaia DR3, we conducted an investigation into the substructures in the stellar halo to obtain their chemical and kinematic information. We employed the HDBSCAN algorithm in the (L z , E, L xy ) space for clustering and utilized Mahalanobis distance to merge the clustering results. Ultimately, within the LAMOST DR9 MRS catalog, we identified members of Gaia–Sausage–Enceladus (GSE), Helmi streams, and the high-α disk and in situ halo. Additionally, there are three clusters that cannot be correlated with known substructures. Based on their angular momentum (L z ) characteristics, we have tentatively named two of them High-L z -Cluster-1 (HLC-1) and High-L z -Cluster-2 (HLC-2), both exhibiting high angular momentum (L z ) features. The third cluster, owing to its extremely low L z , is termed the Low-L z -Cluster. Then detailed analysis about those substructures are executed including member stars’ identification and kinematic and chemical distribution. The resulting GSE member stars exhibit a well-defined metallicity distribution function with its peak at [Fe/H] ∼ −1.2, consistent with previous findings, and we found there appear to be some differences in the distribution on the V total–r gal plane between GSE and its included Splash component. Regarding Helmi streams, in the [Mg/Fe]–[Fe/H] plane we detected potential traces of the α-knee phenomenon at [Fe/H] ∼ −1.53. Finally, we tentatively provide the possible origin for the HLC-1 and HLC-2.

Exploring the origin of the distance bimodality of stars in the periphery of the Small Magellanic Cloud with APOGEE and Gaia

ABSTRACT The Magellanic Cloud system represents a unique laboratory for study of both interacting dwarf galaxies and the ongoing process of the formation of the Milky Way and its halo. We focus on one aspect of this complex, three-body interaction – the dynamical perturbation of the Small Magellanic Cloud (SMC) by the Large Magellanic Cloud (LMC), and specifically potential tidal effects on the SMC’s eastern side. Using Gaia astrometry and the precise radial velocities (RVs) and multielement chemical abundances from Apache Point Observatory Galactic Evolution Experiment (APOGEE-2) Data Release 17, we explore the well-known distance bimodality on the eastern side of the SMC. Through estimated stellar distances, proper motions, and RVs, we characterize the kinematics of the two populations in the bimodality and compare their properties with those of SMC populations elsewhere. Moreover, while all regions explored by APOGEE seem to show a single chemical enrichment history, the metallicity distribution function (MDF), of the ‘far’ stars on the eastern periphery of the SMC is found to resemble that for the more metal-poor fields of the western periphery, whereas the MDF for the ‘near’ stars on the eastern periphery resembles that for stars in the SMC Centre. The closer eastern periphery stars also show RVs (corrected for SMC rotation and bulk motion) that are, on average, approaching us relative to all other SMC populations sampled. We interpret these trends as evidence that the near stars on the eastern side of the SMC represent material pulled out of the central SMC as part of its tidal interaction with the LMC.

The JWST Resolved Stellar Populations Early Release Science Program. IV. The Star Formation History of the Local Group Galaxy WLM

We present the first star formation history (SFH) and age–metallicity relation (AMR) derived from resolved stellar populations imaged with the JWST NIRCam instrument. The target is the Local Group star-forming galaxy WLM at 970 kpc. The depth of the color–magnitude diagram (CMD) reaches below the oldest main sequence turnoff with a signal-to-noise ratio = 10 at M F090W = + 4.6 mag. This is the deepest CMD for any galaxy that is not a satellite of the Milky Way. We use Hubble Space Telescope (HST) optical imaging that overlaps with the NIRCam observations to directly evaluate the SFHs derived based on data from the two great observatories. The JWST and HST-based SFHs are in excellent agreement. We use the metallicity distribution function measured from stellar spectra to confirm the trends in the AMRs based on the JWST data. Together, these results confirm the efficacy of recovering an SFH and AMR with the NIRCam F090W−F150W filter combination, and validate the sensitivity and accuracy of stellar evolution libraries in the near-infrared relative to the optical for SFH recovery work. From the JWST data, WLM shows an early onset to star formation, followed by an extended pause post-reionization before star formation reignites, which is qualitatively similar to what has been observed in the isolated galaxies Leo A and Aquarius. Quantitatively, 15% of the stellar mass formed in the first Gyr, while only 10% formed over the next ∼5 Gyr. The stellar mass then rapidly doubled in ∼2.5 Gyr, followed by constant star formation over the last ∼5 Gyr.

Metallicity Distribution Functions of 13 Ultra-faint Dwarf Galaxy Candidates from Hubble Space Telescope Narrowband Imaging

We present uniformly measured stellar metallicities of 463 stars in 13 Milky Way (MW) ultra-faint dwarf galaxies (UFDs; M V = −7.1 to −0.8) using narrowband CaHK (F395N) imaging taken with the Hubble Space Telescope. This represents the largest homogeneous set of stellar metallicities in UFDs, increasing the number of metallicities in these 13 galaxies by a factor of 5 and doubling the number of metallicities in all known MW UFDs. We provide the first well-populated MDFs for all galaxies in this sample, with 〈[Fe/H]〉 ranging from −3.0 to −2.0 dex, and σ [Fe/H] ranging from 0.3–0.7 dex. We find a nearly constant [Fe/H]∼ −2.6 over 3 decades in luminosity (∼102–105 L ⊙), suggesting that the mass–metallicity relationship does not hold for such faint systems. We find a larger fraction (24%) of extremely metal-poor ([Fe/H]< −3) stars across our sample compared to the literature (14%), but note that uncertainties in our most metal-poor measurements make this an upper limit. We find 19% of stars in our UFD sample to be metal-rich ([Fe/H] > −2), consistent with the sum of literature spectroscopic studies. MW UFDs are known to be predominantly >13 Gyr old, meaning that all stars in our sample are truly ancient, unlike metal-poor stars in the MW, which have a range of possible ages. Our UFD metallicities are not well matched to known streams in the MW, providing further evidence that known MW substructures are not related to UFDs. We include a catalog of our stars to encourage community follow-up studies, including priority targets for ELT-era observations.

Transiting Exoplanet Yields for the Roman Galactic Bulge Time Domain Survey Predicted from Pixel-level Simulations

The Nancy Grace Roman Space Telescope (Roman) is NASA’s next astrophysics flagship mission, expected to launch in late 2026. As one of Roman’s core community science surveys, the Galactic Bulge Time Domain Survey (GBTDS) will collect photometric and astrometric data for over 100 million stars in the Galactic bulge in order to search for microlensing planets. To assess the potential with which Roman can detect exoplanets via transit, we developed and conducted pixel-level simulations of transiting planets in the GBTDS. From these simulations, we predict that Roman will find between ∼60,000 and ∼200,000 transiting planets—over an order of magnitude more planets than are currently known. While the majority of these planets will be giants (R p > 4R ⊕) on close-in orbits (a < 0.3 au), the yield also includes between ∼7000 and ∼12,000 small planets (R p < 4R ⊕). The yield for small planets depends sensitively on the observing cadence and season duration, with variations on the order of ∼10%–20% for modest changes in either parameter, but is generally insensitive to the trade between surveyed area and cadence given constant slew/settle times. These predictions depend sensitively on the Milky Way’s metallicity distribution function, highlighting an opportunity to significantly advance our understanding of exoplanet demographics, in particular across stellar populations and Galactic environments.

Exploring the Sun’s birth radius and the distribution of planet building blocks in the Milky Way galaxy: a multizone Galactic chemical evolution approach

ABSTRACT We explore the influence of the Milky Way galaxy’s chemical evolution on the formation, structure, and habitability of the Solar system. Using a multizone Galactic chemical evolution (GCE) model, we successfully reproduce key observational constraints, including the age–metallicity ([Fe/H]) relation, metallicity distribution functions, abundance gradients, and [X/Fe] ratio trends for critical elements involved in planetary mineralogy, including C, O, Mg, and Si. Our GCE model suggests that the Sun formed in the inner Galactic disc, Rbirth,⊙ ≈ 5 kpc. We also combined a stoichiometric model with the GCE model to examine the temporal evolution and spatial distribution of planet building blocks within the Milky Way galaxy, revealing trends in the condensed mass fraction (fcond), iron-to-silicon mass fraction (firon), and water mass fraction (fwater) over time and towards the inner Galactic disc regions. Specifically, our model predicts a higher fcond in the protoplanetary disc within the inner regions of the Milky Way galaxy, as well as an increased firon and a decreased fwater in the inner regions. Based on these findings, we discuss the potential impact of the Sun’s birth location on the overall structure and habitability of the Solar system.

The stellar halo in Local Group Hestia simulations

Stellar chemical abundances and kinematics provide key information for recovering the assembly history of galaxies. In this work we explore the chemo-chrono-kinematics of accreted and in situ stellar populations, by analyzing six M31/Milky Way (MW) analogues from the HESTIA suite of cosmological hydrodynamics zoom-in simulations of the Local Group. We show that elemental abundances ([Fe/H], [Mg/Fe]) of merger debris in the stellar haloes are chemically distinct from the survived dwarf galaxies, in that they are [α/Fe]-enhanced and have lower metallicity in the same stellar mass range. Therefore, mergers debris have abundances expected for stars originating from dwarfs that had their star formation activity quenched at early times. Accreted stellar haloes, including individual debris, reveal [Fe/H] and [Mg/Fe] gradients in the E − Lz plane, with the most metal-rich, [α/Fe]-poor stars, which have formed in the inner parts of the disrupted systems before the merger, contributing mainly to the central regions of the host galaxies. This results in negative metallicity gradients in the accreted components of stellar haloes at z = 0, seen also for the individual merger debris. We suggest, therefore, that abundance measurements of halo stars in the inner MW will allow constraining better the parameters, such as the mass and merger time, of MW’s most massive merger Gaia-Sausage-Enceladus. The metallicity distribution functions (MDFs) of the individual debris show several peaks and the majority of debris have lower metallicity than the in situ stars in the prograde part of the E − Lz space. At the same time, non-rotating and retrograde accreted populations are very similar to the in situ stars in terms of [Fe/H] abundance. Prograde accreted stars show a prominent knee in the [Fe/H]–[Mg/Fe] plane, reaching up to solar [Mg/Fe], while retrograde stars typically contribute to the high-[Mg/Fe] sequence only. We find that the most metal-poor stars ([Fe/H] ≲ −1) of the HESTIA galaxies exhibit net rotation up to 80 km s−1, which is consistent with the Aurora population recently identified in the MW. At higher metallicities ([Fe/H] ≈ −0.5 ± 0.1) we detect a sharp transition (spin-up) from the turbulent phase to a regular disk-like rotation. Different merger debris appear similar in the [Fe/H]–[Mg/Fe] plane, thus making it difficult to identify individual events. However, combining a set of abundances, and especially stellar age, makes it possible to distinguish between different debris.

The energy distribution of the first supernovae

ABSTRACT The nature of the first Pop III stars is still a mystery and the energy distribution of the first supernovae is completely unexplored. For the first time we account simultaneously for the unknown initial mass function (IMF), stellar mixing, and energy distribution function (EDF) of Pop III stars in the context of a cosmological model for the formation of a MW-analogue. Our data-calibrated semi-analytic model is based on a N-body simulation and follows the formation and evolution of both Pop III and Pop II/I stars in their proper time-scales. We discover degeneracies between the adopted Pop III unknowns, in the predicted metallicity and carbonicity distribution functions and the fraction of C-enhanced stars. None the less, we are able to provide the first available constraints on the EDF, $dN/dE_\star \propto E_{\star }^{-\alpha _e}$ with 1 ≤ αe ≤ 2.5. In addition, the characteristic mass of the Pop III IMF should be mch &amp;lt; 100 M⊙, assuming a mass range consistent with hydrodynamical simulations (0.1–1000 M⊙). Independent of the assumed Pop III properties, we find that all $\rm [C/Fe]\gt +0.7$ stars (with $\rm [Fe/H]\lt -2.8$) have been enriched by Pop III supernovae at a $\gt 20~{{\ \rm per\ cent}}$ level, and all $\rm [C/Fe]\gt +2$ stars at a $\gt 95~{{\ \rm per\ cent}}$ level. All very metal-poor stars with $\rm [C/Fe]\lt 0$ are predicted to be predominantly enriched by Pop III hypernovae and/or pair instability supernovae. To better constrain the primordial EDF, it is absolutely crucial to have a complete and accurate determination of the metallicity distribution function, and the properties of C-enhanced metal-poor stars (frequency and [C/Fe]) in the Galactic halo.

Modeling the Chemical Enrichment History of the Bulge Fossil Fragment Terzan 5

Terzan 5 is a heavily obscured stellar system located in the inner Galaxy. It has been postulated to be a stellar relic, a bulge fossil fragment witnessing the complex history of the assembly of the Milky Way bulge. In this paper, we follow the chemical enrichment of a set of putative progenitors of Terzan 5 to assess whether the chemical properties of this cluster fit within a formation scenario in which it is the remnant of a primordial building block of the bulge. We can explain the metallicity distribution function and the runs of different element-to-iron abundance ratios as functions of [Fe/H] derived from optical-infrared spectroscopy of giant stars in Terzan 5 by assuming that the cluster experienced two major star formation bursts separated by a long quiescent phase. We further predict that the most metal-rich stars in Terzan 5 are moderately He-enhanced, and we predict a large spread of He abundances in the cluster, Y ≃ 0.26–0.335. We conclude that current observations fit within a formation scenario in which Terzan 5 originated from a pristine or slightly metal-enriched gas clump about one order of magnitude more massive than its present-day mass. Losses of gas and stars played a major role in shaping Terzan 5 the way we see it now. The iron content of the youngest stellar population is better explained if the white dwarfs that give rise to type Ia supernovae (the main Fe factories) sink toward the cluster center, rather than being stripped by the strong tidal forces exerted by the Milky Way in the outer regions.

Light-curve Recovery with the Vera Rubin Observatory’s LSST. I. Pulsating Stars in Local Group Dwarf Galaxies

This investigation is mainly focused on the LSST Survey Strategy Optimization process, a bottom-up approach that turned out to be quite effective in involving the scientific community in the definition of the LSST observing strategy. We are mainly interested in using radial variables (RR Lyrae, classical Cepheids, long-period variables) as stellar tracers and distance indicators, and we developed a new tool called PulsationStarRecovery to quantify the recovery of the light-curve period and amplitude from an LSST-simulated time series. The outputs of this code are pulsation parameters (period, amplitude, mean magnitude) together with quantitative information concerning the difference between the shape of the light curve and template light curves. Furthermore, we apply the newborn metric to simulate LSST observations and recovery of different types of pulsating stars hosted by selected massive stellar systems (19 Local Group dwarf galaxies and the Large Magellanic Cloud) to show how the recovery changes according to distance and variable-star type. We show that this exercise is essential to understand the potential of LSST in this field since excellent recovery is necessary to optimize the use of predicted period–luminosity, period–amplitude, and color–color relations to constrain the cosmic distance scale and the metallicity distribution function of different stellar populations.

Detailed α abundance trends in the inner Galactic bulge

Context. Until now, heavy interstellar extinction has meant that only a few studies of chemical abundances have been possible in the inner Galactic bulge. However, it is crucial to learn more about this structure in order to better understand the formation and evolution of the centre of the Galaxy and galaxies in general. Aims. In this paper, we aim to derive high-precision α-element abundances using CRIRES high-resolution IR spectra of 72 cool M giants of the inner Galactic bulge. Methods. Silicon, magnesium, and calcium abundances were determined by fitting a synthetic spectrum for each star. We also incorporated recent theoretical data into our spectroscopic analysis (i.e. updated K-band line list, better broadening parameter estimation, non-local thermodynamic equilibrium (NLTE) corrections). We compare these inner bulge α abundance trends with those of solar neighbourhood stars observed with IGRINS using the same line list and analysis technique; we also compare our sample to APOGEE DR17 abundances for inner bulge stars. We investigate bulge membership using spectro-photometric distances and orbital simulations. We construct a chemical-evolution model that fits our metallicity distribution function (MDF) and our α-element trends. Results. Among our 72 stars, we find four that are not bulge members. [Si/Fe] and [Mg/Fe] versus [Fe/H] trends show a typical thick disc α-element behaviour, except that we do not see any plateau at supersolar metallicities as seen in other works. The NLTE analysis lowers [Mg/Fe] typically by ∼0.1 dex, resulting in a noticeably lower trend of [Mg/Fe] versus [Fe/H]. The derived [Ca/Fe] versus [Fe/H] trend has a larger scatter than those for Si and Mg, but is in excellent agreement with local thin and thick disc trends. With our updated analysis, we constructed one of the most detailed studies of the α abundance trends of cool M giants in the inner Galactic bulge. We modelled these abundances by adopting a two-infall chemical-evolution model with two distinct gas-infall episodes with timescales of 0.4 Gyr and 2 Gyr, respectively. Conclusions. Based on a very meticulous spectral analysis, we have constructed detailed and precise chemical abundances of Mg, Si, and Ca for cool M giants. The present study can be used as a benchmark for future spectroscopic surveys.

Unravelling the mass spectrum of destroyed dwarf galaxies with the metallicity distribution function

ABSTRACTAccreted stellar populations are comprised of the remnants of destroyed galaxies, and often dominate the ‘stellar haloes’ of galaxies such as the Milky Way (MW). This ensemble of external contributors is a key indicator of the past assembly history of a galaxy. We introduce a novel statistical method that uses the unbinned metallicity distribution function (MDF) of a stellar population to estimate the mass spectrum of its progenitors. Our model makes use of the well-known mass–metallicity relation of galaxies and assumes Gaussian MDF distributions for individual progenitors: the overall MDF is thus a mixture of MDFs from smaller galaxies. We apply the method to the stellar halo of the MW, as well as the classical MW satellite galaxies. The stellar components of the satellite galaxies have relatively small sample sizes, but we do not find any evidence for accreted populations with L &amp;gt; Lhost/100. We find that the MW stellar halo has N ∼ 1−3 massive progenitors (L ≳ 108L⊙) within 10 kpc, and likely several hundred progenitors in total. We also test our method on simulations of MW-mass haloes, and find that our method is able to recover the true accreted population within a factor of 2. Future data sets will provide MDFs with orders of magnitude more stars, and this method could be a powerful technique to quantify the accreted populations down to the ultra-faint dwarf mass scale for both the MW and its satellites.

A Photometric Survey of Globular Cluster Systems in Brightest Cluster Galaxies

Hubble Space Telescope imaging for 26 giant early-type galaxies, all drawn from the MAST archive, is used to carry out photometry of their surrounding globular cluster (GC) systems. Most of these targets are Brightest Cluster Galaxies and their distances range from 24–210 Mpc. The catalogs of photometry, completed with DOLPHOT, are publicly available. The GC color indices are converted to [Fe/H] through a combination of 12 Gyr single stellar population models and direct spectroscopic calibration of the fiducial color index (F475W–F850LP). All of the resulting metallicity distribution functions (MDFs) can be accurately matched by bimodal-Gaussian functions. The GC systems in all of the galaxies also exhibit shallow metallicity gradients with projected galactocentric distance that average . Several parameters of the MDFs including the means, dispersions, and blue/red fractions are summarized. Perhaps the most interesting new result is the trend of blue/red GC fraction with galaxy mass, which connects with predictions from recent simulations of GC formation within hierarchical assembly of large galaxies. The observed trend reveals two major transition stages: for low-mass galaxies, the metal-rich (red) GC fraction f(red) increases steadily with galaxy mass, until halo mass M h ≃ 3 × 1012 M ⊙. Above this point, more than half the metal-poor (blue) GCs come from accreted satellites and f(red) starts declining. But above a still higher transition point near M h ≃ 1014 M ⊙, the data hint that f(red) may start to increase again because the metal-rich GCs also become dominated by accreted systems.

Spectroscopic follow-up of statistically selected extremely metal-poor star candidates from GALAH DR3

ABSTRACT The advent of large-scale stellar spectroscopic surveys naturally leads to the implementation of machine learning techniques to isolate, for example, small sub-samples of potentially interesting stars from the full data set. A recent example is the application of the t-SNE statistical method to ∼600 000 stellar spectra from the GALAH survey in order to identify a sample of candidate extremely metal-poor (EMP; [Fe/H] ≤ –3) stars. We report the outcome of low-resolution spectroscopic follow-up of 83 GALAH EMP candidates that lack any previous metallicity estimates. Overall, the statistical selection is found to be efficient (∼one-third of the candidates have [Fe/H] ≤ –2.75) with low contamination (&amp;lt;10 per cent have [Fe/H] &amp;gt; –2), and with a metallicity distribution function that is consistent with previous work. Five stars are found to have [Fe/H] ≤ –3.0, one of which is a main sequence turnoff star. Two other stars are revealed as likely carbon-enhanced metal-poor (CEMP) stars of type CEMP-s, and a known carbon star is re-identified. The results indicate that the statistical selection approach employed was successful, and therefore it can be applied to forthcoming even larger stellar spectroscopic surveys with the expectation of similar positive outcomes.

The multizone chemical evolution of the Galactic bulge: predicting abundances for different radial zones

ABSTRACT Due to its proximity, the stellar populations of the Galactic bulge (GB) can be resolved and can be studied in detail. This allows tracing the bulge metallicity distribution function (MDF) for different spatial regions within the bulge, which may give us clues about the bulge formation and evolution scenarios. In this work, we developed a chemical evolution model (CEM), taking into account the mass distribution in the bulge and disc, to derive the radial dependence of this time-scale in the Galaxy. Since the infall rate depends on that time-scale in the CEM, the results of the model were used to test a scenario where the bulge is formed inside-out. The obtained results for the [α/Fe] versus [Fe/H] relationship, the MDF and the [Fe/H] radial gradient in the bulge have been compared to available data in the literature. The model is able to reproduce most of the observational data: the spread in the relation [α/Fe] versus [Fe/H], the MDF shape in different regions of the bulge, the [Fe/H] radial gradient inside it, and the age–metallicity relation, as well as the [α/Fe] evolution with age. The results of the model point to a scenario where the bulk of the bulge stars pre-existed the boxy/peanut X-shape bar formation. As a result, the classical origin of the GB is not ruled out, and this scenario may be invoked to explain the chemical properties of the GB.

A method for identifying metal-poor stars with Gaia BP/RP spectra

Context. The study of the oldest and most metal-poor stars in our Galaxy promotes our understanding of the Galactic chemical evolution and the beginning of Galaxy and star formation. However, they are notoriously difficult to find, with only seven stars at [Fe/H] &lt; −5.0 having been detected to date. Thus, the spectrophotometric data of 219 million sources which became available in the third Gaia Data Release comprise a very promising dataset for the identification of metal-poor stars. Aims. We want to use the low-resolution Gaia Blue Photometer/Red Photometer (BP/RP) spectra to identify metal-poor stars. Our primary aspiration is to help populate the poorly constrained tail of the metallicity distribution function of the stellar halo of the Galaxy. Methods. We developed a metal-poor candidate selection method based on flux ratios from the BP/RP Gaia spectra, using simulated synthetic spectra. Results. We found a relation between the relative iron abundance and the flux ratio of the Ca H &amp; K region to that of the Hβ line. This relation is temperature and surface gravity dependent, and it holds for stars with 4800 K ≤ Teff ≤ 6300 K. We applied it to noisy simulated synthetic spectra and inferred [Fe/H] with an uncertainty of σ[Fe/H] ≲ 0.65 dex for −3 ≤ [Fe/H] ≤ 0.5 and G = 15–17 mag, which is sufficient to identify stars at [Fe/H] &lt; −2.0 reliably. We predict that by selecting stars with inferred [Fe/H] ≤ −2.5 dex, we can retrieve 80% of the stars with [Fe/H] ≤ −3 and have a success rate of about 50%, that is one in two stars we select would have [Fe/H] ≤ −3. We do not take into account the effect of reddening, so our method should only be applied to stars which are located in regions of low extinction.

TOPoS

Context. The TOPoS project has the goal to find and analyse Turn-Off (TO) stars of extremely low metallicity. To select the targets for spectroscopic follow-up at high spectral resolution, we have relied on low-resolution spectra from the Sloan Digital Sky Survey. Aims. In this paper we use the metallicity estimates we have obtained from our analysis of the SDSS spectra to construct the metallicity distribution function (MDF) of the Milky Way, with special emphasis on its metal-weak tail. The goal is to provide the underlying distribution out of which the TOPoS sample was extracted. Methods. We make use of SDSS photometry, Gaia photometry and distance estimates derived from the Gaia parallaxes to derive a metallicity estimate for a large sample of over 24 million TO stars. This sample is used to derive the metallicity bias of the sample for which SDSS spectra are available. Results. We determined that the spectroscopic sample is strongly biased in favour of metal-poor stars, as intended. A comparison with the unbiased photometric sample allows to correct for the selection bias. We select a sub-sample of stars with reliable parallaxes for which we combine the SDSS radial velocities with Gaia proper motions and parallaxes to compute actions and orbital parameters in the Galactic potential. This allows us to characterize the stars dynamically, and in particular to select a sub-sample that belongs to the Gaia-Sausage-Enceladus (GSE) accretion event. We are thus able to provide also the MDF of GSE. Conclusions. The metal-weak tail derived in our study is very similar to that derived in the H3 survey and in the Hamburg/ESO Survey. This allows us to average the three MDFs and provide an error bar for each metallicity bin. Inasmuch the GSE structure is representative of the progenitor galaxy that collided with the Milky Way, that galaxy appears to be strongly deficient in metal-poor stars compared to the Milky Way, suggesting that the metal-weak tail of the latter has been largely formed by accretion of low mass galaxies rather than massive galaxies, such as the GSE progenitor.

Probing the Galactic halo with RR lyrae stars − III. The chemical and kinematic properties of the stellar halo

ABSTRACT Based on a large spectroscopic sample of ∼4300 RR Lyrae stars with metallicity, systemic radial velocity, and distance measurements, we present a detailed analysis of the chemical and kinematic properties of the Galactic halo. Using this sample, the metallicity distribution function (MDF) as a function of r and the velocity anisotropy parameter β profiles (for different metallicity populations) are derived for the stellar halo. Both the chemical and kinematic results suggest that the Galactic halo is composed of two distinct parts, the inner halo and outer halo. The cutoff radius (∼30 kpc) is similar to the previous break radius found in the density distribution of the stellar halo. We find that the inner part is dominated by a metal-rich population with extremely radial anisotropy (β ∼ 0.9). These features are in accordance with those of ‘Gaia-Enceladus-Sausage’ (GES) and we attribute this inner halo component as being dominantly composed of stars deposited from this ancient merged satellite. We find that GES probably has a slightly negative metallicity gradient. The metal-poor populations in the inner halo are characterized as a long tail in MDF with an anisotropy of β ∼ 0.5, which is similar to that of the outer part. The MDF for the outer halo is very broad with several weak peaks and the value of β is around 0.5 for all metallicities.